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Abstract 

We consider a large number of two-level atoms interacting with the mode of a cavity in the 
rotating- wave approximation (Tavis-Cummings model). We apply the Holstein-Primakoff trans- 
formation to study the model in the limit of the number of two-level atoms, all in their ground 
state, becoming very large. The unitary evolution that we obtain in this approximation is applied 
to a macroscopic superposition state showing that, when the coherent states forming the superpo- 
sition are enough distant, then the state collapses on a single coherent state describing a classical 
radiation mode. This appear as a true dynamical effect that could be observed in experiments with 
cavities. 
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I. INTRODUCTION 



The foundations of quantum optics rely on the Hamiltonian of the interaction of a single 
radiation mode with an atom that, for all practical applications in optical regime, reduces 
to the well-known Jaynes-Cummings model Q]. The proper working of this model is due 
to the effectiveness of certain approximations [2j: The rotating wave approximation and 
the two-level approximation. These approximations are enough to justify the success of the 
Jaynes-Cummings model in the optical regime. 

In experiments involving such a model a quite common effect that appears is decoherence, 
ihat is the loss of unitarity in the quantum evolution due to external environmental effects 



3, 



Decoherence can make very difficult to realize quantum computers. Indeed, our aim in 
this paper is to treat the problem of quantum coherence for a system obeying the Jaynes- 
Cummings model and being composed by N two-level atoms in the limit of N becoming 
infinitely large. Such a model is known in the literature as the Tavis-Cummings model 

Recent works on decoherence have enlarged the meaning of this mechanism to a general 
reduction of the density matrix to the mixed form removing interference without any need 
of coupling to an external environment. The possible existence of such kind of decoherence 
without environment is, indeed, an open question. 

Decoherence can also be seen as an intrinsic effect due to the unitary evolution when 
the „ Um be r of particles beco m es ve ry , arg e fl fl g A firS t e^ple in this sense has 
been given by Gea-Banacloche Q that computed asymptotic states for the radiation field 
in the Jaynes-Cummings model in the limit of a large photon number. Such states have 
been recently observed in experiments by Haroche et al. 10, ^| proving the emergence of 
classical behavior already with very few photons. The work of Gea-Banacloche appears in 
some way dual to the one we present here. 

In order to study the Tavis-Cummings Hamiltonian in the limit of a very large number 
of particles on a macroscopic superposition state we will adopt an approach firstly devised 
by Persico and Vetri and recently applied in a paper by Berman et al. ^| . These 
authors use the Holstein-Primakoff 1^ transformation to change the contribution of the 
two-level atoms into a bosonic field. In this way, one gets the Hamiltonian of two coupled 
harmonic oscillators that can be easily diagonalized. The main result of this paper will be 
to show how the Tavis-Cummings Hamiltonian can produce decoherence, intrinsically in 
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this approximation, making this effect observable in experiments with cavities. The unitary 



evolution is applied to a superposition state of two coherent states of radiation. W 



ren the 
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distance between such states is large enough, as already seen in recent experiments 
decoherence sets in and the superposition disappears leaving just a coherent state. 

The paper is structured in the following way. In seclH] we introduce the model and 
the Holstein-Primakoff transformation. In sec lllll the higher order corrections to the unitary 
evolution are evaluated. In sec II VI the unitary evolution of a macroscopic superposition state 
is obtained proving that decoherence is produced when the components of the superposition 
state are significantly distant. In secElthe conclusions are given. 



II. TAVIS-CUMMINGS MODEL AND HOLSTEIN-PRIMAKOFF TRANSFOR- 
MATION 

The N-atom Jaynes-Cummings model was firstly considered by Cummings and Tavis 4] 
and can be written as 

^ N N 

H = ua^a + — ^2a a i + g^2(<T +i a + a- i a^) (1) 
1 i=i i=i 

being uj the frequency of the radiation mode, a and a' the ladder operators, A the separation 
between the energy levels of the two level atoms, g the coupling and a^,a + i and a_j the 
Pauli spin matrices. By introducing the operators S x = | Y<iLi S± = Z)i=i °±i we can 
rewrite the above Hamiltomian in the form 

H = ua)a + AS Z + g{S + a + S-cfi). (2) 

A "ferromagnetic state" is characterized by having all the two-level atoms in their ground 
state. For the sake of simplicity we put us at the resonance A = uj. Our aim is to study the 
Tavis-Cummings model in this situation. This is accomplished by the so called Holstein- 
Primakoff transformation. This is defined by 

S + = &t(_2S-&t 6 )* ( 3 ) 

S- = (-25-6 f 6) 3 6 
S z = S + b ] b 
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being S = — y and bosonic operators. The aim of this transformation is to obtain 



a series in the parameter being |i?| 3> 1, that is pertinent to our situation. These 



computations are well-known in literature 



but we report it here for completeness: 



H = + uo{a)a + tfb) + VNg{(Jb + tfa) (4) 

Hi = ?=(a j b j bb + atftfb) 



and the second and third order corrections are given by 

H 2 = --^={a!tfbtfbb + atftfbtfb) (5) 

H 3 = ^^(a^tfbtfbtfbb + atftfbtfbtfb) 

48v / A^ 

from which an easy rule to obtain higher order terms is realized. So, n-th order term is given 
by the multiplicative factor — q n — ^-r, being q n the numerical coefficient of the corresponding 
order in the series of (1 — x)a, and Hamilton operator is given by a) TYk^ftb^b + h.c. 
Our approximation holds until the average number of bosonic excitations described by the 
operator b is largely smaller than \S\. 

The leading order Hamiltonian can be immediately diagonalized by introducing two 
bosonic operators c\ and c 2 as 

a — b 

that give 

H = (u + VNgylc! + {u- ^Ng)c\c 2 (7) 

i.e. two independent harmonic oscillators and we have omitted the constant. This means 
that, at the leading order, we expect the coherence to be preserved In sec JIVI we will 
show that this conclusion can be evaded in some way. Anyhow, the unitary evolution at the 
leading order is given by 

Uo(t) = exp — it(u + vNg)c\ci exp —it(u — vNg)G 2 C2 . (8) 

The eigenvalues are given by 

e nin2 (u,g,N) = n 1 (uj + VNg) + n 2 (uj-\^Ng). (9) 



The eigenstates are given by 



\nl-n2) = ^=r|0). 



>n 2 \ 



■|o>- 



(10) 



having put |0)_ the state having all the two- level atoms in the ground state on which the c\ 
operator is acting. This clarifies how the time evolution happens in the limit of N two-level 
atoms all in their ground state. The next order correction is O 
For a coherent state of the radiation mode 



|^(0)) = e QaWa |0) 



N 

T 



(ii) 



we will get at the leading order 



\^{t)) = |5e- <(wW ^ )t ) + |ae- <(w -^ )t ). 



(12) 



with a = being 



c\\ae 



)- 



ae 



i{LO+VNg)t I ^ e -i{co+VNg)t 



>- 



(13) 



and 



C2 | 5e -i("-^)*)_ = fie -<(«-^)t| 5e -i( W -V3v fl )t^_ (M) 
coherent states of the two harmonic oscillators. These equations give the useful result for 



our aims 



a\^(t)) = ae~ luJt cos(V Ngt)\ilj(t)). 
From this equation it is easy to recover the well-known results 
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(n) = (4){t)\a ] 'a\ip{t)) = \a\ 2 cos 2 (VNgt) 



being h the number operator for the radiation field and 



(15) 



(16) 



(n 2 ) - (n) 2 = \a\ 2 cos 2 (VNgt) 



(17) 



that proves that at the leading order the coherence is kept. This results holds, at this order, 
also in the thermodynamic limit iV — > oo, g — > and yNg = constant. 
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III. HIGHER ORDER CORRECTIONS TO THE UNITARY EVOLUTION 



The answer to the question about the coherence being kept to higher order has been 
properly answered in Ref . [13I] . The answer is given by computing the first order correction 
to the eigenvalues by the Rayleigh-Schrodinger equation. The Hamiltonian to use is given 
by Hi that is, using the operators c\ and c 2 

Hl = - l(clci)2 " (4c2) !| 4c2 " clci " [cW ~ 4c2)(c!c2 + Cl4) + c{a - ~ 4cJ (18) 

and we obtain the correction [13[ 
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4VN 



[n\ -n 2 2 + n 2 - m) 



that gives the following correction to the average of the number operator 



(h) 



\a\ 



, ,9 \n>\ 2ni \n\ 2n2 
' J2 e- 2 '"' 1 1 , 11 , cos 



Ngt + 
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WN 



>i + n 2 )t) 



(19) 



(20) 



ri\\ n 2 \ 

and we lose coherence if the second term in the argument of the cosine is comparable to the 
first one but, in the thermodynamic limit as defined above, we are granted that coherence 
is maintained. 

It is interesting to compute also the correction to the eigenstates (fTUj). Rayleigh- 
Schrodinger series gives 



K; n 2 )i = — (nix/mV^ + IK - 1; n 2 + 1) + n 2 ^ 2 ^Jn^+\\n x + 1; n 2 - 1)) (21) 
oiV 

that is O (jjj while, as seen above, the correction to the eigenvalues is O • Then, we 
are able to compute the correction to the average of the number operator to first order. This 
gives 



6(h) 



1 ^ QlAia |a| 2 "i Ul 2 ^ 



E ■ 

ni,n 2 



-2|a| 2 



nil 712! 



{(m + n 2 ){nl + ml) + n\{n 2 + 1) + n\[n x + 1) (22) 



+ {[(m - l) 2 + (n 2 + l) 2 ]^^n 2 + 1 + [(m + l) 2 + (n 2 - l) 2 ]\/ni + 1^} x 



cos 



2VNgt + 
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2VN 



(ni + n 2 - l)t) 



+ \(n 2 + 2)y/TnVn 2 + lVni - Wn 2 + 2 + (m + 2)^^ + 1V^2 - V^i + 2] x 



cos 



WNgt + -j=(ni + n 2 - l)t) 



and again we have a confirmation that the thermodynamic limit grants that coherence is kept 
for an initial coherent state. As our aim is to apply the time evolution to a superposition of 
coherent states, this conclusion is crucial for our result to hold in the thermodynamic limit. 



IV. UNITARY EVOLUTION OF A MACROSCOPIC SUPERPOSITION STATE 

Now, let us consider a phase Schrodinger cat state [j| 

\M0)) = ^{\ f y^*) + \^~ 4 *)) (23) 

being 

1 1 

A/ 2 = - - (24) 

2 1 + cos( 7 2 sin(20))e- A2 ' K ' 

A 2 = 2j 2 sin 2 the distance between the coherent states, 7 and tf> two real numbers. The 
time evolution of this state is given by 

\tfj s (t)} = Af r|^e _<[(w+ ^ fl) *~* 1 ) + |7e~ i[(w_v ^ fl)t_ ^)_ + \^ e ~ i[ ^ + ^ 9)t+ ^) + \^e~ i[ ^~^ 9)t+ ^)- 

(25) 

with 7 = The computation of the averages is rather straightforward and gives using 
eq.flHJ) 

, , , 9 9/ /T7 , 1 + cos(20 + 7 2 sin(20))e~ A . . 

{Mt)WMMt)) = i^A^N'A T^^hm^ (26) 

that has the property to give the result of a single coherent state or (f> — > but, mostly 
important, for A ^> 1 one recover the same result of a single coherent state. So, more distant 
are the coherent states and easier the superposition is removed in the time evolution. This 
is confirmed by the computation: 

toWlaWohM*)) = 7 4 cosHVNgt) 1 + ™ s W + 7* ^ ^ (27) 



+ 7 2 cos 2 (v / iV^t) 



1 + cos( 7 2 sin(20))e A 
1 + cos(20 + 7 2 sin(20))e- A 



1 + cos(7 2 sin(20))e A ' 2 

that gives back the result of a single coherent state both for = and A ^> 1 confirming that 
the unitary evolution for the Tavis-Cummings model, in the thermodynamic limit, grants 
the disappearance of the superposition state. This situation has been encountered also for 
the Dicke model in Ref.Q] where it was shown an identical result in this case but resorting 
to the concept of singular limits. Here, the argument still relies on the thermodynamic limit, 
that is essential to the proof, but we have also assumed the components of the Schrodinger 
Celt StcltjG are enough distant, A>1. 

It is important to note that, if the coherent states in the superposition are truly macro- 
scopic, we are left with a single coherent state behaving as a classical radiation field. 
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V. CONCLUSIONS 



In conclusion, we have shown how decoherence can be produced in the thermodynamic 
limit also for the N atom Jaynes-Cummings model that can be straightforwardly used to 
test the very existence of a somewhat different kind of decoherence, even if the N atoms can 
be seen as an environment. But decoherence appears dynamically. 

We would like to emphasize the duality of our approach with respect to the one of 
Gea-Banacloche [$]. Striking evidence of classicality emerging by increasing the number 
of photons in a cavity, without resorting to an y c oncept of environment, has been given 
recently in the experiment of Haroche's group giving full support to the analysis of 
Gea-Banacloche |9]. So, it would be interesting to see also states for N atoms in the ther- 
modynamic limit as obtained here. Then, decoherence can also appear as an intrinsic effect 
of the unitary evolution in the thermodynamic limit. 
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